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The idea of sensitivity in ant colony systems has been exploited in hybrid ant-based models with promising results 
for many combinatorial optimization problems. Heterogeneity is induced in the ant population by endowing individual 
ants with a certain level of sensitivity to the pheromone trail. The variable pheromone sensitivity within the same 
population of ants can potentially intensify the search while in the same time inducing diversity for the exploration of 
the environment. The performance of sensitive ant models is investigated for solving the generalized vehicle routing 
problem. Numerical results and comparisons are discussed and analysed with a focus on emphasizing any particular 
aspects and potential benefits related to hybrid ant-based models. 
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1. Introduction 

The potential of ant-based models [1, _I |_} in solving 
difficult optimization problems has been well emphasized 
by successful results obtained in many and varied fields 
including transportation optimization, quadratic assign- 
ment, scheduling, vehicle routing and protein folding. In- 
spired by the real-world collective behaviour of social in- 
jects, Ant Colony System (ACS) algorithms [2j rely on 
the stigmergic interactions between many identical artifi- 
cial ants to find solutions to a given problem. Each ant 
generates a complete tour (associated to a problem solu- 
tion) by probabilistically choosing the next node at each 
path intersection based on the cost and the amount of 
pheromone on the connecting edge. Stronger pheromone 
trails are preferred and the most promising tours build up 
higher amounts of pheromone in time. 

Inducing heterogeneity in the population by enabling 
each artificial ant to react in a different way to the same 
environment [loj represents a promising approach to the 
application of ant-based models for solving complex real- 
world problems possibly with a dynamic character. Each 
individual ant can be endowed with a certain level of sensi- 
tivity to the pheromone trail triggering various types of re- 
actions to a changing environment. The variable pheromone 
sensitivity within the same population of ants can poten- 
tially intensify the search (normally through high sensi- 
tivity levels) while in the same time inducing diversity for 
the exploration of the environment. The decision of a low- 
level sensitive ant regarding the action to be performed 



Email addresses: cmpintea@cs.ubbcluj.ro (Camelia-M. 
Pintea), chira@cs.ubbcluj.ro (Camelia Chira), 
ddumitrScs.ubbcluj.ro (D. Dumitrcscu), pop_petricaSyahoo.com 
(Petrica C. Pop) 



crucially contributes to the quality of the search process 
and solutions. 

The Generalized Vehicle Routing Problem (GVRP) is 
an extension of the Vehicle Routing Problem (VRP) and 
was introduced by Ghiani and Improta @. The GVRP 
is the problem of designing optimal delivery or collection 
routes, subject to capacity restrictions, from a given depot 
to a number of predefined, mutually exclusive and exhaus- 
tive node-sets (clusters). 

The GVRP belongs to the class of generalized combina- 
torial optimization problems, which are natural extensions 
of combinatorial optimization problems by considering a 
related problem relative to a given partition of the nodes 
of the graph into node sets, while the feasibility constraints 
are expressed in terms of the clusters. In the literature we 
can find several generalized problems such as the general- 
ized minimum spanning tree problem (see [l_|), the gen- 
eralized traveling salesman problem, the generalized vehi- 
cle routing problem, the generalized (subset) assignment 
problem, etc. These generalized problems belong to the 
class of NP-completc problems, are harder than the clas- 
sical ones and nowadays are intensively studied due to the 
interesting properties and applications in the real world, 
even though many practitioners are reluctant to use them 
for practical modeling problems because of the complexity 
of finding optimal or near-optimal solutions. 

Ghiani and Improta @ showed that the problem can 
be transformed into a capacitated arc routing problem 
(CARP) and Baldacci et al. proved that the reverse 
transformation is valid. Recently, Pop (l2j provided a new 
efficient transformation of the GVRP into the classical ve- 
hicle routing problem (VRP). As far as we know, the only 
specific algorithm for solving the GVRP was developed by 
Pop et al. 11| and was based on ant colony optimization. 



The aim of this paper is to investigate the performance 
of the Sensitive Ant Model (SAM) [10| in solving the Gen- 
eralized Vehicle Routing Problem (GVRP) We report 
numerical results of the SAM model for several GVRP 
benchmark problems and discuss the performance of SAM 
compared to the standard ACS technique. 

2. Definition and Complexity of the GVRP 

Let G = (V,A) be a directed graph with V = {0, 1,2, 
as the set of vertices and the set of arcs A — \ i,j € 

V,i 7^ j}. A nonnegative cost Cij associated with each arc 
(i, j) £ A. The set of vertices (nodes) is partitioned into 
fc+1 mutually exclusive nonempty subsets, called clusters, 
V ,Vi, V k (i.e. V = V U Vt U ... U V k and Vj n V p = for 
all {0, 1, k} and I ^ p). The cluster Vo has only 

one vertex 0, which represents the depot, and remaining 
n nodes belonging to the remaining k clusters represent 
geographically dispersed customers. Each customer has a 
certain amount of demand and the total demand of each 
cluster can be satisfied via any of its nodes. There exist 
m identical vehicles, each with a capacity Q. 

The generalized vehicle routing problem (GVRP) con- 
sists in finding the minimum total cost tours of starting 
and ending at the depot, such that each cluster should be 
visited exactly once, the entering and leaving nodes of each 
cluster is the same and the sum of all the demands of any 
tour (route) does not exceed the capacity of the vehicle Q. 
An illustrative scheme of the GVRP and a feasible tour is 
shown in the Figure [T] 



m=2 and Q=25 



q,=12 




Figure 1: An example of a feasible solution of the GVRP 

The GVRP reduces to the classical Vehicle Routing 
Problem (VRP) when all the clusters are singletons and 
to the Generalized Traveling Salesman Problem (GTSP) 
when m — 1 and Q = oo. 



The GVRP is AP-hard because it includes the gener- 
alized traveling salesman problem as a special case when 
m = 1 and Q = oo. 

Several real- world situations can be modelled as a GVRP. 
The post-box collection problem described in Laportc ct 
al. [7| becomes an asymmetric GVRP if more than one ve- 
hicle is required. Furthermore, the GVRP is able to model 
the distribution of goods by sea to a number of customers 
situated in an archipelago as in Philippines, New Zeeland, 
n } Indonesia, Italy, Greece and Croatia. In this application, 
a number of potential harbours is selected for every island 
and a fleet of ships is required to visit exactly one harbour 
for every island. 

Several applications of the GTSP (Laporte et al. [|j]) 
may be extended naturally to GVRP. In addition, several 
other situations can be modeled as a GVRP, these include: 

• the Traveling Salesman Problem (TSP) with profits 
(Feillet et al. 0); 

• a number of Vehicle Routing Problem (VRP) exten- 
sions: the VRP with selective back hauls, the cover- 
ing VRP, the periodic VRP, the capacitated general 
windy routing problem, etc.; 

• the design of tandem configurations for automated 
guided vehicles (Baldacci et al. [l[). 

3. The ACS algorithm for solving GVRP 

The ACS-based algorithm for GVRP [H| uses artificial 
ants in order to construct vehicle routes by successively 
choosing exactly one node from each cluster. This task 
continues until each cluster has been visited. Whenever 
the choice of another node from a cluster would lead to an 
infeasible solution because of vehicles capacity, the depot 
is chosen and a new route is started. 

Initially, the Nearest Neighbor (NN) algorithm - with 
the rule always go to the nearest as-yet-unvisited location - 
is considered. The best solution of Nearest Neighbor (NN) 
algorithm (L + ) is used for ACS-based algorithm start. 

The number of ants corresponds to the number of GVRP 
customers m. At the beginning of an iteration, an ant is 
placed at each node (customer). After initializing the ba- 
sic ant system algorithm, the two steps: (i) construction 
of vehicle routes and (ii) trail update are repeated for a 
given number of iterations. 

To favor the selection of an edge with a high pheromone 
level and high visibility, a probability function pi is de- 
fined as follows: 
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where is the set of unvisited neighbors of node i by ant 
k, j S J,- and 8 is a parameter used for tuning the relative 
importance of visibility. 
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After an artificial ant has constructed a feasible solu- 
tion, the pheromone trails are laid depending on the ob- 
jective value Lfc. For each edge that was used by ant k, 
the pheromone trail is updated according to the following 
rule: 



Table 1: Problem characteristics for the ant-based algorithms for 
GVRP 
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Tij (t + 1) = (l-p) nj (t)+ pjj: 



(2) 



where p 6 (0, 1) is an evaporation rate parameter. 

A tabu list prevents ants visiting clusters they have 
previously visited. The ant tabu list is cleared after each 
completed tour. 

The global update rule, applied by the elitist ants, as 
in ACS H is: 



i (t+l) = (l-p)n j (t)+p 



1 



(3) 



where L + is the so far best solution. 

4. Sensitive Ant-based Model for GVRP 

The Sensitive Ant Model (SAM) technique proposed 
in is engaged for solving the GVRP. The general ap- 
proach to solving GVRP using SAM is the same with the 
ACS approach presented in the previous section except 
that the transition probabilities defined by SAM are used. 
The initialization of the algorithm, the update rules and 
the maintenance of a tabu list are kept the same in SAM 
for GVRP. 

The SAM algorithm involves several ants able to com- 
municate in a stigmergic manner (influenced by pheromone 
trails) for solving complex search problems. Within the 
SAM model, each ant is characterized by a pheromone sen- 
sitivity level (PSL). The PSL value is expressed by a real 
number in the unit interval [0, 1]. When PSL is null the 
ant completely ignores stigmergic information and when 
PSL is one the agent has maximum pheromone sensitiv- 
ity. The ants with a low PSL value are more independent 
and are considered environment explorers. They have the 
potential to autonomously discover new promising regions 
of the solution space. The ants with high PSL values are 
very sensitive to pheromone traces. They are influenced 
by stigmergic information and therefore intensively exploit 
the promising search regions already identified. 

SAM introduces a measure of randomness proportional 
to the level of individual PSL in the decisions of ants re- 
garding the path to follow. This is achieved by modifying 
the transition probabilities using the PSL values in a renor- 
malization process [10] . The SAM renormalized transition 
probability for ant k (influenced by PSL) is denoted by 
spl (t) and is given by the following equation: 



sp»(t)=p»(t)-PSL k (t), 



(4) 



where p*£ (t) is the probability for ant k to choose the 
next node j from current node i (as given in ACS - see 
Equation [1]) and PSL)~(t) represents the PSL value of ant 
k at time t. 
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llcil51 


2 


160 


320 
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16eil76A 
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140 


280 


10 
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16cil76B 
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100 


300 


15 
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16cil76C 
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180 


360 
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16eil76D 
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220 


440 
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21eill01A 


2 


200 


400 
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21eill01B 


2 


112 


224 


14 


7 



It can be noticed that the SAM probability of select- 
ing the next node is the same with the ACS one when 
PSL value is one. In order to associate a standard prob- 
ability distribution to the system, the SAM virtual state 
corresponding to the 'lost' probability of (1 — PSLk(t)) 
has to be defined. The associated virtual state decision 
rule specifies the action to be taken when the virtual state 
is selected using the renormalized transition mechanism. 
The following rule is used in the current paper: the ant 
randomly chooses an available node with uniform proba- 
bility if the virtual state is selected. This approach favors 
the increasing of randomness in the selection process with 
the decreasing of sensitivity level to pheromone. 

5. Computational results 

The performance of the SAM and ACS for solving 
GVRP is investigated. Numerical experiments focus on 
seven benchmark problems from the TSPLIB library 14 1 . 
These problems contain between 51 and 101 customers 
(nodes) , which are partitioned into a given number of clus- 
ters, and in addition the depot. 

Originally the set of nodes in these problems is not 
divided into clusters. The CLUSTERING procedure pro- 
posed by Fischetti et al. 5] is used to divide data into 
node-sets. This procedure sets the number of clusters 
m = [j], identifies the m farthest nodes from each other 
and assigns each remaining node to its nearest center. 

Table[T]contains the description of the GVRP instances 
addressed in this paper. 

The meaning associated with the columns in Table [T] is 
as follows: 

• Problem: The name of the test problem contains the 
number of clusters (first digits in the problem name) 
and the number of nodes (last digits in the problem 
name). 

• VR: The minimal number of vehicles needed for a 
route in order to cover even the largest capacity of a 
cluster (VR= Vehicles/Route) 

• Q': the capacity Q ■ VR, where Q is the capacity of 
a vehicle available at the depot. 
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Table 2: Best Values and Times - ACS and SAM algorithms for 
solving GVRP 



Problem 


ACS 


Time ACS 


SAM 


Time SAM 


llcilSl 


418.85 


212 


418.21 


297 


16eil76A 


668.78 


18 


651.98 


25 


16eil76B 


625.83 


64 


599.23 


166 


16cil76C 


553.21 


215.00 


577.49 


88 


16cil76D 


508.81 


177.00 


515.64 


120 


21eill0lA 


634.74 


72 


634.74 


111 


21eill01B 


875.58 


8.00 


966.17 


52 



Table 3: Average Values and Times - ACS and SAM algorithms for 
solving GVRP 



Problem 


ACS 


Time ACS 


SAM 


Time SAM 


lleil51 


429.85 


210.20 


424.05 


96.40 


16eil76A 


706.09 


109.20 


677.85 


187.60 


16eil76B 


684.04 


50.7 


608.62 


173.7 


16cil76C 


625.87 


73.00 


602.06 


42.25 


16eil76D 


566.56 


93.20 


533.12 


223.45 


21eill01A 


699.46 


29.00 


690.39 


124.30 


21eill01B 


996.41 


25.95 


998.71 


27.90 



The same parameter setting was used in both SAM 
and ACS algorithms in order to allow a meaningful direct 
comparison: tq — 0.1 (the initial value of all pheromone 
trails), a = l,P = 5,p= 0.0001 and q = 0.5. In the SAM 
algorithm, the PSL value is randomly generated between 
and 1 for each ant. 

Numerical results indicate a competitive performance 
of the SAM algorithm. Figure[5]presents SAM results from 
20 successive runs for the considered problem instances. 



SAM Results 
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Figure 2: SAM results from 20 runs for the considered problem in- 
stances Il4l 

Tables [2] and [3] present comparative numerical results 
obtained in 20 runs. The performance of SAM in solving 



GVRP is compared to that of ACS [11 



The following information is contained in Tables[5]and|3J 

• Best length: the minimal length of collection routes; 

• Best time: the time of the minimal collection routes; 

• Avg. length: the average length for 20 runs; 

• Avg. time: the average time (in seconds) for 20 runs. 

The computational values are the result of the aver- 
age of 20 successively runs of both algorithms. Termina- 
tion criterion is either the maximum number of iterations, 
Nuer = 250000 or the maximum running time (five min- 
utes) on a AMD 2600, 1.9Ghz and 1024 MB. 



When comparing the best solution reported in 20 runs, 
the performance of the ACS and SAM algorithms are sim- 
ilar. SAM reports better values for three out seven prob- 
lems while ACS does better for other three problems. It 
should be noticed that whenever a method is able to ob- 
tain a better solution, it also reports a longer running time 
compared to the other one. 

The average solutions obtained by SAM are clearly im- 
proved compared to ACS although the running time has 
slightly increased for some of the GVRP instances. SAM 
detects a better average solution (calculated based on the 
25 runs) for six out of the seven benchmark problems. 
Overall better average SAM results are facilitated by a 
better exploration of the search space and exploitation of 
new solutions. This is due to the variable sensitivity in- 
duced in SAM via random individual PSL values. 

6. Conclusions 

Sensitive heterogeneous ant-based models facilitate a 
balanced search process by endowing ants with different 
pheromone sensitivity levels translated into different search 
strategies. An effective exploration of the search space is 
performed particularly by ants having low pheromone sen- 
sitivity while the exploitation of intermediary solutions is 
facilitated by highly-sensitive ants. 

The performance of hybrid ant-based models is investi- 
gated with successful results for solving the NP-hard Gen- 
eralized Vehicle Routing Problem. Variable pheromone 
sensitivity in ant-based models proves to be benefic to the 
search process leading to better results compared to the 
ant colony system algorithm. Numerical results encour- 
age the exploration of new ways to induce heterogeneity 
in ant-based models. 
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